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We demonstrate the feasibility of computing a 3 from the lattice three-gluon vertex in the Landau gauge. Data 
from 16 4 and 24 4 quenched lattices at (3 — 6.0 are presented. Our main result is that 2-loop asymptotic scaling 
is observed for momenta in the range 1.8 — 2.3 GeV, where lattice artifacts appear to be under control. 



1. THE METHOD 

a s can be extracted from the 3-gluon vertex jlj 
by evaluating two-point and three-point off-shell 
gluon Green's functions on the lattice and impos- 
ing non-perturbative renormalisation conditions 
on them, for different values of the external mo- 
menta. The main advantages of this method are 
the possibility to obtain a s (/z) at several momen- 
tum scales n from a single simulation and the fact 
that lattice perturbation theory (LPTH) is not 
needed to match our coupling to a—. Also, the 
method can be applied with no modifications to 
the unquenched case. On a more technical level, 
additional advantages are that, since one works 
in momentum space, lattice artifacts can be care- 
fully analysed, as shown is section 2.1, and since 
only gluonic operators are used, the technology 
is simple. The method is fully described in 
Working in the Landau gauge and using a suit- 
able definition of the lattice gluon field A^, one 
can define the momentum space gluon propaga- 
tor G { J\jp) = T^{jp) Gu(p 2 ), where T^(p) is 
the transverse projector, and the complete gluon 
three-point function G\ J ' al3l (p 1 ,p 2 ,P3)- The 
propagator is non-perturbatively renormalised by 
imposing that for p 2 — /i 2 it attains its continuum 



tree-level value: 

G R (p)\ p 2 =ll 2 = Z A 1 {^a)G u {pa)\ p 2 = 



(1) 



The three-gluon vertex is evaluated at the 
asymmetric kinematical points where, based on 
the general form of continuum the vertex func- 
tion in the Landau gauge |J, the following rela- 
tion holds: 

,4 ^(3) 

1:111 ! I "-1)1!: 

§iZy X (pa) g pp. (2) 



(Gu(pa)) 2 Gu(0) 

By computing the above ratio one can determine 
the factor Zy 1 go, where Zy is the vertex renor- 
malisation constant. Finally, one can define the 
running coupling g at the scale [i from the renor- 
malised three-gluon vertex as 

5(m) = Z 3 A /2 (tm) Zy x (iia) g . (3) 

and a s (fj.) = g(fi) 2 /4n. This choice corresponds 
to a momentum^subtraction scheme, usually re- 
ferred to as MOM in continuum QCD M. 

2. NUMERICAL RESULTS 

We consider two quenched data sets at j3 = 6.0; 
150 configurations on a 16 4 lattice, and 103 con- 
figurations on a 24 4 lattice. We used a Landau 
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gauge fixing overrelaxation algorithm, with the 
final iterations being performed in double preci- 
sion. It is worth emphasising that the numerical 
accuracy of the gauge-fixing is crucial to obtain a 
good signal for the three-point function. 



2.1. Analysis of Lattice Artifacts 

A comparison of numerical values of ratios of 
Green's functions tensor components with the 
theoretical expectations from LPTH and con- 
tinuum perturbation theory (CPTH) provides a 
check for discretisation errors. In the case of the 
gluon propagator, the numerical values are in per- 
fect agreement with LPTH. As for the three-gluon 
vertex, in Table 1 we show various ratios of ten- 
sor components on the 16 4 lattice, where in the 
first column the lattice momentum components 
are specified. Unless otherwise noted, the uncer- 
tainty is less than one unit in the last quoted fig- 
ure. Only the ratio Gioi/Gm is in poor agree- 
ment with LPTH. This is related to the fact that 
the Landau gauge condition does not fix this ra- 
tio. Overall, the tensor structure is the one ex- 
pected from LPTH. This enables us to estimate 
the size of the lattice artifacts and to identify a 
"continuum window" in momentum space. 

2.2. Renormalisation Constants and Run- 
ning Coupling 

In Figure |we plot <?(/z) on the 16 4 lattice. All 
data are plotted vs. [i = expressed in GeV. 

In order to detect violations of rotational invari- 
ance, we have used whenever possible different 
combinations of lattice momentum vectors for a 
fixed value of p 2 , plotting separately the corre- 
sponding data points. 

We obtain a clear signal, although the data 
show some violation of rotational invariance. 

The important question is whether there is a 
momentum range where our coupling runs ac- 
cording to the two-loop expression 
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= boH 



A 2 ' 

MOM 



h , , , y 2 
V lnln{ A*7 
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where b a = 11/16tt 2 , bi = 102/(16tt 2 ) 2 and 
A.-^tt, is the QCD scale parameter for our 

MOM 

scheme. To investigate this point, we compute 
A„-rrr, as a function of the measured values of 

MOM 



Figure 1. Running coupling g(^) vs. /x on the 
smaller lattice. 



g 2 (/i) according to the formula 



A m^m = » ex P 
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2b g 2 (fi) 



[b g 2 ^)} ^.(5) 



If the coupling runs according to 



then 



^mom as defined above approaches a constant 
value when /i — » oo. By plotting ^-^^f M versus \i 
(see Figure 2), one notices that 




Figure 2. A,-~, vs. u for our two lattice sizes. 

b MOM 



for /j, < 1.8 GeV, ^-^Qi, depends strongly on 
fi. However, in the range 1.8 < ji < 2.3 GeV 
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Table 1 

Symmetry tests for G^- (p, 0, -p) on the 16 4 lattice at f3 = 6.0 (C=CPTH, L=LPTH, N=numerical). 
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the data are consistent with a constant value for 
A „-^rr , . No violations of rotational invariance are 

MOM 

observed and a comparison of the two lattice sizes 
shows no volume dependence. For p > 2.3 GeV, 
rotational invariance is broken by higher order 
terms in a 2 and the two-loop behaviour disap- 
pears. In summary, we appear to have a "contin- 
uum window" in the range 1.8 < p < 2.3 GeV, 
where two-loop scaling is observed and lattice ar- 
tifacts are under control. In order to extract a 
prediction for A^^g-^, we fit the data points in 
the continuum window to the curve obtained by 
inserting (^) in @ . We take as our best estimate 
the fit to the 16 4 data, for which the statistical 
errors are smaller, and obtain 



A M^M=°- 



I ± 0.02 ± 0.09 GeV, 



(6) 



where the first error is statistical and the second 
error comes from the uncertainty on the value of 
a" 1 . 

3. MATCHING TO MS 

We can extract a prediction for a-jj^ with zero 
quark flavours from our numerical results for 
A.-rx-;,. The ratio Axr^/A,-—, can be deter- 

MOM . MSI MOM 

mined to all orders in the coupling constant from 
a one- loop continuum calculation. We obtain, 
in the Landau gauge and for zero quark flavours 
h^—/h^X~, = 0.35. This result is in agreement 

MS' MOM b 

with previous one-loop calculations of the three- 
gluon vertex ||. Using (JsJ) and the above result, 
we obtain 



A^L = 0.31 ± 0.05 GeV. 

MS 



(7) 



This is the main result of our computation, which 
is in very good agreement with the one ofref. §. 



In terms of aj^L, our result yields: 
aH.(2.0 GeV) = 0.24 ± 0.02. 



(8) 



Although LPTH is not needed for matching, 
we can use it to perform some interesting cross- 
checks of the continuum calculation (see for 
details). 

4. CONCLUSIONS 

We have shown in the quenched approxima- 
tion that a non-perturbative determination of the 
QCD running coupling can be obtained from first 
principles by a lattice study of the three-gluon 
vertex. We have some evidence that systematic 
lattice effects are under control in our calcula- 
tion. LPTH is not needed to match our results 
to MS and the extension to the full theory does 
not present in principle any extra problem. 

C. Parrinello and DGR acknowledge support 
from PPARC through travel grant GR/L29927 
and Advanced Fellowships. We also acknowl- 
edge support from EPSRC grants GR/K41663 
and GR/K55745. 

REFERENCES 

1. C. Parrinello, Phys. Rev. D50 (1994) 4247. 

2. B. Alles et al., Nucl. Phys. B502 (1997) 325. 

3. J.S. Ball, TW. Chiu, Phys. Rev. D22 (1980) 
2550. 

4. A. Hasenfratz, P. Hasenfratz, Phys. Lett. B93 
(1980) 165. 

5. F.T. Brandt and J. Frenkel, Phys. Rev. D33 
(1986) 464. 

6. G.S. Bali and K. Schilling, Phys. Rev. D47 
(1993) 661. 



